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All-optical feedback can be effected by putting the output of a source cavity through a Faraday 
isolator and into a second cavity which is coupled to the source cavity by a nonlinear crystal. If the 
driven cavity is heavily damped, then it can be adiabatically eliminated and a master equation or 
quantum Langevin equation derived for the first cavity alone. This is done for an input bath in an 
arbitrary state, and for an arbitrary nonlinear coupling. If the intercavity coupling involves only the 
intensity (or one quadrature) of the driven cavity, then the effect on the source cavity is identical 
to that which can be obtained from electro-optical feedback using direct (or homodyne) detection. 
If the coupling involves both quadratures, this equivalence no longer holds, and a coupling linear in 
the source amplitude can produce a nonclassical state in the source cavity. The analogous electro- 
optic scheme using heterodyne detection introduces extra noise which prevents the production of 
^ . nonclassical light. Unlike the electro-optic case, the all-optical feedback loop has an output beam 

' (reflected from the second cavity). We show that this may be squeezed, even if the source cavity 

remains in a classical state. 

D . PACS numbers: 42.50.Lc, 42.50.Dv, 03.65.Bz 
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I. INTRODUCTION 



Electro-optic feedback is the use of a photocurrent to control the source of the light incident on the detector 
producing that current. Such feedback has long been used to control noise in optical systems such as lasers. In many 
cases, the noise which is being controlled is classical noise. That is, noise which is well above the shot-noise limit (also 
' known as the quantum limit). However, technological advance now enables experimentalists to work with light which 
is at, or even (thanks to squeezing p|) below, the quantum limit. Thus, the quantum limits to electro-optic feedback 
have been of considerable interest in the past decade 0, 0> 0, H| • Our interest lies with feedback onto the dynamics 
of the source cavity 0, 0, • We have described such feedback using quantum measurement theory, in the form of 
stochastic quantum trajectories 0, 0, 0, 0. In the limit that the time delay in the feedback loop is much less 
Oh! that the cavity lifetime, a master equation describing the feedback can be derived. This is a simple and elegant way 
■ to treat feedback. There is also a corresponding Langevin equation approach 8], which is more convenient for some 
purposes. The clearest result of our theory is that, unless the system without feedback has nonclassical dynamics, then 
controlling its dynamics via a photocurrent cannot produce nonclassical light. That is, feedback based on external 
photodetection cannot produce squeezing. 

Given these limitations on the noise-reduction abilities of electro-optic feedback, it is natural to ask, what about 
all-optical feedback. That is, instead of detecting the emitted light, it is reflected around a loop back into the source 
cavity, or, more fruitfully, into another cavity which is coupled to the first cavity in some way. Of course, for this 
mechanism to be considered feedback, the feedback loop must be one-way, otherwise we would simply be describing 
a pair of doubly- coupled cavities. The general configuration under consideration is shown schematically in Fig. 1. 
One mechanism for achieving the required unidirectionality is the Faraday isolator which utilizes Faraday rotation 
and polarization-sensitive beam splitters. A quantum theoretical treatment which incorporates this spatial symmetry 
breaking at the level of the Hamiltonian has been given recently by Gardiner |13| . If the propagation time between the 
source system the and driven system is ignored, then a master equation for both systems may be derived. This result 
was obtained simultaneously by Carmichael [lj > who introduced the term "cascaded systems" . The new feature which 
we consider here is for the driven system to be coupled back to the source system via some interaction Hamiltonian, 
giving rise to all-optical feedback. 

The quantum theory of cascaded systems is summarized in Sec. II, and the simplest example of feedback via this 
method is considered. In Sec. Ill, we consider a more complex feedback model utilizing a nonlinear crystal to influence 
the source cavity, with the strength of the feedback being proportional to the intensity of the source cavity. Both 
a master equation and a quantum Langevin equation are derived, and are shown to be approximately the same as 
those pertaining to electro-optic feedback via direct detection. Similarly, the quadrature-sensitive all-optical feedback 



FIG. 1: [See end.] Diagram of the general all-optical feedback scheme considered. The annihilation operators for the source and 
driven cavities are denoted ci and C2 respectively, while 62 and 63 represent traveling waves. The nonlinear coupling between 
the cavities is indicated by V, FR denotes a Faraday rotator, and PBS a polarization-sensitive beam splitter. 
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scheme considered in Sec. IV is equivalent to homodyne-mcdiated feedback in its effect on the source cavity. It is 
thus not possible to produce a nonclassical state in the source cavity via a coupling which is linear in the source 
cavity amplitude or intensity (which correspond to classical driving or detuning). Unlike the electro-optic case, the 
all-optical feedback loop has an output beam of light (reflected off the driven cavity), and we show that this may 
be squeezed even though the source cavity is not. Sec. V treats an all-optical feedback scheme which is sensitive to 
both quadratures of the source cavity. This has no direct electro-optic equivalent, and may produce a nonclassical 
source state even with a coupling linear in the source amplitude. An electro-optic analogue can be defined, and the 
extra noise which it introduces seen explicitly. The two schemes are contrasted using the simplest feedback example 
introduced in Sec. II. 



II. CASCADED SYSTEMS 
A. Unidirectional coupling 



In this secti on, we present a summary of the theory of cascaded open quantum systems. For different presentations, 
see Refs. 0, Q. Our starting point is the input-output theory of open quantum systems 0,0]. This theory 
describes a system interacting locally with a bath consisting of a continuum of harmonic oscillators. Physically, the 
system may be an optical cavity, and the bath the external electromagnetic field modes with momentum aligned 
to the cavity axis. The electric field (or rather, one polarization component) at a particular point in space-time 
(paramaterized by z,t) is represented approximately by the Heisenberg picture operator [l6| 

E(z,t) = ^^L[b(z,t) + b\z,t)]. (1) 

Here, A is the cross sectional area of the beam, and only frequencies near the central wave number k are assumed to 
be of interest. The canonical commutation relations for the complex amplitudes b(z,t) are 

[b(z,t),tf(z',t)] = c5(z-z% (2) 

where c is the speed of light, and for regions where the field propagates freely, 

b(z,t + r) = b(z- cr,t). (3) 

Let the external field be coupled to the cavity by a very good mirror at z = 0. The field with z < then represents 
an incoming field and that with z > an outgoing one. We assume a linear coupling of the form 

H x {t) = iVnP> + (0,*)ci(t) - c\(t)b(0,t)}, (4) 

where 71 is the coupling constant having the dimension of inverse time, and C\ (t) is the annihilation operator of the 
cavity tuned to the frequency ck. Ignoring other dynamics, the evolution of an arbitrary Heisenberg operator a(t) is 

d(t) = -[tf(Q,t)a(t) - c\(t)b(0,i),a(t)}. (5) 

Now because of the singularity of the canonical commutation relations J5J , it is necessary to be careful in dealing with 
this evolution equation. A convenient method is to use quantum Ito stochastic differential calculus |l5l Il7| . Define 
an input field, representing the field just before it interacts with the cavity at time t by 

6i(t)=6(0-,t). (6) 

This can be thought of as a white noise term, independent of the state of the cavity at time t. The analog of the 
Weiner increment in the Ito calculus is then 

dBi(i) =bi(t)dt, (7) 

which satisfies 

[dB l {t),dB\{t)]= dt. (8) 
The evolution of an arbitrary operator is then given by 

a{t + dt) = U\(t,t + dt)a(t)U!(t,t + dt) (9) 
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where 



Ux{t,t + dt) =exp{v7T dB\{t)ci{t) - dB x (t)ci{t) } 



(10) 



In Eq. El the bath operators dB\(t) and dB\(t) are independent of the system operator a(t), which makes it an easy 
equation to solve. However, the price which must be paid by this simplification of Ito calculus is that U\{t,t + dt) 
must be expanded to second order in the increment dBi(t). Now if bi(t) is to be thought of as a bath, it should 
be specifiable simply by its moments. We need only the first and second order moments. The first order moment 
corresponds to a coherent amplitude 



(dBi(t)) = f3(t)dt, 



while the second order moments indicate white noise: 



{dB\{t)dB x {i)) = Ndt = (dBi(t)dB{(t)) - 1 
(dB^dB^t)) = Mdt = {dB\(t)dB^ '(*))* . 



(11) 



(12) 
(13) 



These equations include the cases of thermal noise (M = 0), and perfectly squeezed white noise [\M\ 2 = N(N + 1)]. 
It is nonclassical only if \M\ > N. Note that even if N — 0, it is necessary to expand Ui(t, t + dt) to second order 
because then dBi{t)dB\{t) — dt, which could be regarded as vacuum noise. The full result (with time arguments 
omitted for convenience) is the quantum Langevin equation 



da 



7i 
2 



(N + l)(2c\ac 1 - ac\c x - c\c x a) + N(2aac\ 



ac\c 



Cic\a) 



+M[c\, [c{,a}} + M*[c u [ci,a]] \ dt - ^y[[dB\ci - dB lC \, 



(14) 



The final, stochastic term in this equation is essential to preserve canonical commutation relations. However, the 
stochastic terms can be ignored when changing from the Heisenberg to the Schrodinger picture and deriving the 
evolution of the density operator for the cavity mode alone. This is found from the relation 



(da(t)) = Tr[dp(t)a], 



(15) 



where the picture (Schrodinger or Heisenberg) is specified by the placement of the time argument. The resulting 
master equation is 



Pit) 



7i 



(N+l)V[d]p + NV[c\]p + — [c\ , [c\ ,p]] + \ [ci , [a , p]} 



+V7l[/3*(*)ci - (3(t)c\,p), 
where we have defined a superoperator T> taking an arbitrary operator as its argument by 



2?[a]p = apa} 



1 t 1 t 



(16) 



(17) 



It might be thought that this master equation is the only product of the above theory which is of any significance, 
since it generates the evolution of the cavity mode. However, experimentally, it is often the light which leaves the 
cavity through the output mirror, rather than the internal state of the cavity, which is of interest. This is specifically 
so for this paper, in which the output of the cavity will be used in the feedback loop. Thus, we need an expression 
for the field leaving the cavity. From Eq. J3J), this is evidently given by 



b 2 (t) = b(0 + ,t) = Ui(t,t + dtfbi(t)Ui(t,t + dt). 



(18) 



To lowest order in dt, this is 



b 2 (t) = b 1 (t) + ^ 1 - Cl (t). 



(19) 



Just as bi(t) is independent of, and so commutes with an arbitrary system operator a(t') at an earlier time t' < t, the 
output field commutes with all system operators at a later time. 
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Now let this field be the input into another cavity with annihilation operator c 2 (t). If the damping rate for this 
cavity is 72, then the Hamiltonian coupling is 



H 2 (t) = ih^2~[tf(cT,t)c 2 (t) - c\{i)b{ar,t% 



(20) 



where cr is the path length between the two cavities. Proceeding as above, the unitary evolution generated by this 
Hamiltonian is 



U 2 (t,t + dt) =exp{V7^ dBl(t - r)c 2 (t) - dB x (t - r)4(*) }, 
where dB2(t) — b2(t)dt. An arbitrary operator in the source or driven cavity obeys the equation 

a(t + dt) = Ul(t,t + dt)Ul{t,t + dt)a(t)U 2 (t,t + dt)U 1 {t,t + dt). 



(21) 



(22) 



Note that U 2 (t,t + dt) commutes with Ui(t,t + dt) because of the finite r, so the ordering in the above equation is 
not significant. Expanding the terms as above gives 



da = y 



(N + l)(2c\aa - ac\cx - c\aa) + N(2ciac\ 



acic\ 



c\c\a) 



72 
" 2 



+ M[c\, [c\,a]} +M*[ci, [c x ,a]] dt 



(N + l)(2c|ac2 — ac2C2 — 4 c 2 a ) + N(2c 2 ac 2 



ac 2 c 2 



c 2 c\a) 



M[4,[cl,a]]+M*[c 2 ,[c 2 ,a]] dt 
y[[dB\ci - dB\c\,a\ - ^i\dB\c 2 - dB 2 c\, a] 



(23) 

Here, the implicit time argument of a, c\ and C2 is t, while that of dB2 = dB\ + y/jic±dt is t — t. 

Equation ll'-il) is an Ito equation in that it is gives an explicit algorithm for calculating an infinitesimal increment 
in some operator, given all of the other operators at the start of the time interval. However, it is not an Ito equation 
in the sense that the noise terms are not independent of the other operators. Although the noise input dB\{t) is 
independent, dB\{t — r) in dB2(t — r) is not independent of an arbitrary system operator a(t), as can be seen as 
follows 



a(t) = a(t - t) - y/^[dB\(t - r) Cl (t - r) - dBx(t-r)c\(t - r),a(t - r)] + O(r). 



(24) 



To derive a master equation, it is necessary to take the formal limit r — > 0. The basic physics of the problem is 
independent of r because, as yet, we are not considering feedback. This means that corrections of order r in Eq. Ij24(l 
can be ignored, after calculating the second order Ito corrections. By using this equation carefully, Eq. (|23|l can be 
converted to an Ito equation in which the noise terms are independent: 



da = (N + 1) — (2c[aci — ac\c\ — c[cia) + — (2c 2 ac 2 — ac\c 2 — c 2 c 2 a) 



c[c 2 a) 



+M 



ac 2 c 2 



+\/7i72(c2 ac i - ac\ci + c\ac 2 

7 1 /'O t t t 1 1 T 2 /n 1 

+ v / 7i72(c2ac{ — ac 2 c{ + c\ac) 2 — cic^a) 
-±[c\, [c\,a]] + y[4, [4, a]] + V7i72[4> [c\, a}} 
— [ci, [ci,a]] + y[c 2 , [c 2 ,a]] + V7i72[ci, [c 2 ,a] 



c 2 c\a) 



7i[c?Sjci - dB lC {, a] - ^[dB\c 2 - dB lC f 2 ,a\. (25) 

In this equation, all operators have the same time argument. It agrees with the results of Refs. (0,^3) f° r ^ ne case 
N = M = 0, but it should be noted that the details of the derivation in both of these papers differs from ours. 

It is now a simple matter to convert this stochastic Heisenberg equation into a master equation for the density 
operator of both systems W 

W 



(N + 1) 7l X>[ci]W + 7 2^[c 2 ]V^ + V7I72" [ciW, 4] + [c 2 , Wc\ 



5 



y [4, [4, w]] + 1[4, [4. w]] + [4, w]i 



+N 
+M 
+M* 



^[ Cl ,[c u W]] + ^[c 2 ,[c 2 ,W\ 



V7i72[c2, [ci, W]] 
^(t)ca-/3(t)4w]-i[ff,W] 



(26) 

Here we have finally included the intrinsic evolution for the two systems, generated by the Hamiltonian HH. This is 
the general equation for two open quantum systems, linked unidirectionally by a bath of harmonic oscillators with an 
optical frequency coherent amplitude contaminated by white noise. (In the following sections, the coherent amplitude 
will be ignored.) It is a generalization of the equations published by Gardiner [13j and Carmichael 14], which treated 
only a bath in the vacuum state. It has the necessary property that, if H is the sum of Hamiltonians operating in 
the Hilbert subspaces of the two systems, then the source system (ci) is unaffected by the driven system {c%). That 
is to say, the density operator p for the source system (obtained by tracing over the driven system) obeys the original 
master equation This is evident from the fact that the only terms in Eq. I|2()|) containing operators from both 

systems involve an exterior commutator with a driven system operator, which gives zero when traced over the driven 
subspace. It is not possible in general to derive a master equation for the second system alone. Its evolution is literally 
driven by the source system. 



B. Feedback 



The simplest imaginable all-optical feedback scheme would be to take the light from one end of a cavity and reflect 
it back into the other, using a Faraday isolator to prevent interference from reflections in the opposite directions. 
For simplicity, we assume a bath in the vacuum state and a cavity with equal transmittivities at both end-mirrors. 
Denote the input vacuum state by &i(t) = v(f), the cavity mode annihilation operator by a, and the damping rate by 
7. Then the evolution of the cavity mode due to the first mirror is 

a{t) = -la{t)-^v{t), (27) 

and the output field is ^/ja + v(t). If the time delay in the loop is r, then the effect of the feedback is 

M*) = -|o(*) - V7[\/7"a(* - r) + u(t - t)}. (28) 

Provided that jt -C 1, the effect of the time delay can be ignored, apart from introducing an arbitrary phase factor 
depending on the optical path length of the loop. Because the bath is in a vacuum state, it is not necessary to worry 
about the Ito corrections used in going from Eq. (|23[l to Eq. (|25H . Thus the total evolution for a is 

a{t) = -(l + e^)[ 1 a{t) + ^ V {t)\. (29) 

The master equation equivalent to this quantum Langevin equation is 

p = 27(1 + cos 4>)V[a]p — sin (fi[a'a, p] . (30) 

Evidently, if <f> = it, the damping through the mirrors can be completely eliminated. Of course, this is an approximation 
only. In reality, the lifetime of the cavity is enhanced by a factor of order I/7T ^> 1. Losses in the loop would also 
decrease the effectiveness of the feedback. For = 0, the damping rate is doubled, and for other values of 0, the 
cavity becomes detuned, as well as having its damping rate altered. In Sec. V we will show that these features can 
be reproduced by feedback using a heterodyne detection photocurrent to control a driving field at the second mirror. 
However, electro-optic feedback will necessarily introduce extra noise. The all-optical model here obviously does not 
introduce noise; its results are a consequence of classical geometrical optics. 



III. INTENSITY FEEDBACK 



In this section we consider feedback of the form described in the introduction. The light from the source cavity 
is used to drive a second cavity, which coincides spatially at least in part with the source cavity (see Fig. 1). The 
two modes interact via a nonlinear crystal, allowing the driven cavity to control the source cavity. The form of the 
interaction is assumed to depend only on the intensity of the driven cavity. It is thus analogous to feedback of the 
photon flux output of the source cavity. We analyze the feedback system using both master and quantum Langevin 
equations. 
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A. Master equation 



The form of our feedback master equation is that of the general cascaded systems equation (1261 derived in the 
preceding section. What distinguishes it as feedback is that the Hamiltonian TiH is assumed to consist of a source 
cavity Hamiltonian hH a plus an interaction between the source and driven cavity of the form 

V = hc\c 2 K (31) 

where K is an Hermitian operator on the source cavity. In this, and all following sections, we will measure time 
in inverse units of the decay constant 71 for the source cavity. The decay constant 72 for the driven cavity will be 
denoted 7. Because of the interaction term, the dynamics of the source cavity is no longer independent of that of the 
driven cavity. Thus we must consider the master equation for the density operator of both modes W . According to 
Eq. i|26|) . this obeys the master equation 

W = -i[H + 4c 2 K, W] + V[ci]W + ~fV[c 2 ]W + yfi ([ Cl W, c\] + [c 2 , Wc\] \ . (32) 

Note that here we have assumed that the bath is in the vacuum state (N = 0). The reason for this is will become 
apparent later. 

In deriving Eq. I|32[l . it is of course necessary to set r = 0. However, this is no longer a formal mathematical limit. 
Rather, it is a physical condition that the time delay in the feedback loop should be negligible. This is usually desirable 
in feedback loops. A substantial time delay may lead to instabilities and chaos, rather than control. Similarly, for the 
feedback to be effective, the second cavity should respond much faster than the first. This ensures that the state of c 2 
is effectively slaved to that of c\ , and the interaction effects instantaneous feedback as far as the source is concerned. 
Thus, in the limit 7 3> 1, it should be possible to derive a master equation including feedback for the source density 
operator p alone. To do this we note that since 7> 1 and the bath is in the vacuum state, the driven cavity will be 
very close to being in the vacuum state also. This enables the adiabatic elimination procedure we have used previously 
[l8| to be employed again. 

We expand W in powers of 1/ y/j as 

W = po ® |0> <0| + [pi (8 |1)<0| + H.c] + P2 <g> |1)(1|, (33) 

where the ps exist in the source subspace and the other operators in the driven subspace. There is another term of 
second order in l /V/7 , but it plays no part in the procedure which follows. Substituting the above expansion into the 
master equation (|32[) gives the following coupled equations 

po = IP2 + V7(Pi c i + c iPi) + £oPo (34) 

pi = --7P1 - iKpi - V7[ciPo + 0(1/7)] + Co pi (35) 

fa = ~7P2 - i[K,p 2 ] - V7(piCi +cip{) +£op2, (36) 

Here C^p = T>[ci]p — i[H , p] . In this approximation, the source density operator is p = po + p 2 which evidently obeys 

p=-i[K,p 2 ]+C p. (37) 

To turn this into a master equation, we require an expression for p 2 in terms of p ~ po. It is now obvious why we 
assumed a zero temperature bath. For N finite, p 2 would have a finite size irrespective of 7. Thus the signal due to 
the driving from the source, which is of order 1/7 would be swamped by the noise, and the feedback would not work. 

To obtain p 2 it is first necessary to obtain pi. Since almost all of the probability is in po, it is evident from Eq. P5|l 
that p\ relaxes much rapidly than p ~ p. It is thus permissible to set p\ equal to its steady state value of 

2K\^ -2 



Pl =\l + l — j — ClP . (38) 

The expression on the right hand side will be a well defined operator if we assume that \K\ <C 7 in some sense. This 
assumption will be valid in practice, as single-photon nonlinearities are typically much smaller than damping rates. 
It allows the denominator to be expanded to first order in K/j. Since p\ is now slaved to po, it will evolve at a rate 
much smaller than 7. Thus, from Eq. (|36|l . p 2 will relax to a steady state determined by the slaved value of py. 



Pi = ~1P2 



i[K,p 2 }+ 4e lP c\- 4ij~ 1 [K,c 1 pc\}+ C p 2 . (39) 
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The slaved value of P2, again to first order in K/j, is 

4cipct 4i 
P2 = 7T 

Substituting this into Eq. 1)37(1 gives the master equation 



AK 



, cipc[ 



where we have defined 



p = -i[Z,cxpc\] - -[Z, [Z,cipc\]] + T>[c 1 ]p- i[H ,p], 



(40) 



(41) 



Z = iK/j. 



(42) 



This master equation is the general equation for Markovian, intensity-dependent feedback in the small Z limit. 
Unfortunately, it is not a valid master equation [Tgj] in the sense that it cannot be written in the form 



p = -i[H,p] + ^2v[c li ]p, 



(43) 



where the are arbitrary operators. However, there is an equation of this form which is equal to Eq. I(41|l when 
expanded to second order in Z . That equation is 

i>= -i[H ,p]+V[e- iZ Cl }p. (44) 

This equation has been previously derived as a feedback master equation appropriate to electro-optical feedback 
This can be seen more clearly by rewriting Eq. ((44(1 as 



pit + dt) = exp 



-iHq - -cjci j dt 



p(t) exp 



iH - ^ c i c i^ dt 



+ e- iZ cxp(t)c\e^dt. 



(45) 



The two terms in this equation can be given an interpretation in terms of density operators conditioned on possible 
photon detections. The norm of each term gives the probability for the associated event. The first term represents the 
conditioned density operator when no photon is detected [lj. Note that it is only infinitesimally changed, and that 
the probability for this event is thus very close to unity. The second term represents the evolved density operator if 
a photon is detected at the time t. This does not happen very often, but when it does, the state of the system jumps 
(changes discontinuously) via the application of the operator C\ [lj. The effect of the feedback is to cause some 
finite unitary evolution immediately following the detection, via the operator e~ %z . Of course this does not change 
the norm of this term. The physical interpretation of the operator Z is that it can be derived from a time-dependent 
feedback Hamiltonian 



H fb = HZI(t), (46) 

where I(t) is the photocurrent derived from the detector, measured in units of detections per second. This could be 
produced by a variety of electro-optic devices. 

We thus see that there is a strong correspondence between electro-optical feedback using direct detection and all- 
optical feedback with the intensity dependent coupling ((31(1 . One conclusion which can be drawn from this analogy 
is that, as was previously established for electro-optic feedback y|, the all-optic feedback considered in this section 
cannot produce a nonclassical state in the source cavity if K is a classical operator. By "classical" in this context 
we mean a linear combination of ci,c{, and c\c%. (K is of course Hermitian.) Given that the nonlinear interaction 
Hamiltonian 131(1 is already of second order in the field of the driven (C2) cavity, any coupling with a nonclassical 
K would require at least a five- wave mixing interaction (including an auxiliary pump field). This seems exceedingly 
impractical. Thus we can conclude that intensity-dependent all-optical feedback is not a practical way to produce a 
nonclassical source state. 



B. Langevin equation 



The above derivations for the all-optical and electro-optical intensity feedback were done using a master equation 
for the density operator. It is possible to use the Langevin approach for both types of feedback, and this yields some 
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extra information. We begin with the all-optical feedback. The quantum Langevin equation equivalent to the master 
equation is 

1 4 ft T t ft 

a = —(2c[aci — ac[ci — c\c\a) + — (2c 2 ac 2 — ac 2 c 2 — c 2 c 2 a) 

-[b\ci - bic{, a] ~ V7[4 C 2 ~ &2<4> a ] + *[^o + clc 2 K,a}. (47) 

Here b\ is the input vacuum field, and b 2 = b\ + c\ is the output field from the first cavity. This field which is fed 
back is to be understood to be at a slightly earlier time than all of the other operators in the above equation. For 
this reason, it commutes with all other operators. If b 2 is moved to the rear (far right) of any operator expression, 
the vacuum noise of b 2 will not contribute to any average, as the annihilation operator will act directly on the bath in 
the vacuum state. Similar remarks hold for moving b\ to the front of any expression. Thus, if we always put the bath 
operators in normal order as described here, then the fact that it is at a slightly earlier time can be ignored. Of course, 
this technique only works because the bath is in the vacuum state. In general, the time delay causes the corrections 
derived in Sec. II. In the remainder of this section, we will always put b 2 in normal order without mentioning this 
explicitly. 

We wish to adiabatically eliminate mode C2. This obeys 

c 2 = - lea - V7^2 - iKc 2 . (48) 

Now the relaxation rate 7 of mode C2 cannot be much greater than the bandwidth of the vacuum fluctuations, which 
are assumed infinite. Hence, it is not strictly possible to slave C2 to the vacuum fluctuations in 62- However, as far as 
mode ci is concerned, vacuum fluctuations restricted to a bandwidth of 7 are still effectively white, and so their is no 
harm in pretending that c 2 is slaved to the original vacuum fluctuations. This allows us to write the slaved value of 
c 2 as 



C2 = 



2iK 

7 



-^2 (49) 
V7 



Substituting this into the equation for a source cavity operator a gives 

d = -(2e|aci — ac\ci — c\cia) — v^[ci, a] + [c\, a\v 

+i-bl[K,a]b 2 - \b\[K,[K,a])b 2 + i[H Q ,a}. (50) 

7 7 

Here we have put b\ — v because the equation is only valid when the input is a vacuum. Substituting in the expression 
for 6 2 gives 

a = c\aci — —ac\ci — —c\cia — [c\, a] + [c\, a]u 



+(c\ + i/t) U[Z, a] - l[Z, [Z, a}]\ (c x + v) + i[H { 



Qjffl]j (51) 

where Z is as in Eq. I|42|) . It is easy to verify by tracing over the bath that this equation is equivalent to the master 
equation l|41() . 

Just as master equation (|41|l was an approximation to a strictly valid master equation, so is the Langevin equation 
H5ip. The details of the following are to be found in Ref. ||. The exact Langevin equation can be derived from the 
following Hamiltonian 

H fb = TiZb\b 2 . (52) 

This turns out to be equivalent to the photocurrent-dependent Hamiltonian l|4t)|) . That is, simply by replacing the 
photocurrent with the output photon flux b\b 2 , it is possible to derive a Langevin picture equation equivalent to the 
full master equation l|44|) : 

d = c\ac\ — —ac\ci — —c\cia — [c\, a] + [c\, a]v 

+ (c\ + i/t) (e zZ ae~ zZ - a) (ci + v) + i[H , a}. (53) 
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To second order in Z, this is equivalent to the all-optical Langevin equation derived above. Unlike Eq. Q51[l. however, 
Eq. (|53|l is a completely valid Langevin equation so that if ai and a 2 are two arbitrary operators, then the equation 
of motion it generates for the product operator aia 2 is equal to what would be obtained from the two operators 
separately, using the Ito rules for quantum stochastic calculus 

d(aia 2 ) = (dai)a 2 + ai{da 2 ) + (dai)(da 2 ). (54) 

Thus it is also possible to find a correspondence between all-optical and electro-optical feedback using quantum 
Langevin equations. 

One property that a Langevin equation has which the master equation lacks is that it simply gives an expression 
for the field reflected from the driven cavity (see Fig. 1). Calling this field 63, we have 

b 3 = b 2 + ^/jc 2 . (55) 

From the adiabatic expression for c 2 , we have 

4iK~\ 



b, = - 



1 - 

7 



(56) 



Apart from a change of sign (due to reflection), this is equal to what would be obtained from the idealized Hamiltonian 

h = e~ lZ b 2 (57) 

to first order in Z . The statements which were made above regarding the inability of a classical feedback operator Z 
to produce a nonclassical source state do not necessarily apply to the output operator & 3 . It may exhibit nonclassical 
features even if the source state is classical. We will not pursue the properties of this output field further in this 
section, but we will in the next section where we consider quadrature-dependent feedback. 



IV. QUADRATURE FEEDBACK 
A. Master equation 

We begin this section on quadrature feedback by reproducing the calculation of Sec. Ill A, but with the interaction 
between the source and driven cavity 

V = h(c 2 + 4) J, (58) 

which is linear in the real quadrature of the driven cavity. Here J is an Hermitian operator in the source cavity. In 
practice, conservation of energy would require at least one other auxiliary field which may be treated classically. With 
the input to the source cavity in the vacuum state as before, and with the the damping rate 7 of the second cavity 
much greater than that of the first (set to unity), it is again possible to expand the combined density operator in 
powers of l/-^/7 as in Eq. I|33|l . In this case, the source cavity operators obey 

Po = IP2 + V7(Pi c i + ciPi) ~ i[Jpi - Pi J] + £oPo (59) 

Pi = -\lPi-i[Jpo + 0(1/7)] -y/l[c 1 p + 0(1/7)] + A)Pi (60) 

P2 = -IP2 - i[Jp\ - PiA ~ V7(Pi c i + c iPi) +£oP2- (61) 
The approximate source density operator p = po + p 2 obeys 

p = -i[J, pi + p\] + C oP . (62) 

Evidently, to derive a ME for p in this case it is only necessary for p\ to be slaved to p = po + 0(1/7). From 
Eq. I|60(l . this is obviously true for large 7, with the slaved value 

2 _ 

pi = -[-y/jcip-iJp]. (63) 

7 
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Here we are keeping only the lowest order in I/7, and assuming that J is of order ^Fy. Substituting this expression 
into Eq. Ijfi2|l gives the master equation 



P = -i[Y, ap + pc\] - i[y, [Y, p}} + V[ci]p - i[H , p], 



where we have defined 



Y 



-23 



V7 



which is of order unity. This master equation can be put into the form required of master equations thus: 



p = V[c\ — iY]p — i 



H + ^(c\Y + Y Cl ),p 



(64) 



(65) 



(66) 



This feedback master equation has been previously derived from a model of electro-optic feedback using homodyne 
detection [H0- If the output of the cavity is subject to unit efficiency homodyne detection of the real quadrature, 
then the state of the cavity conditioned on the measured photocurrent obeys 



d Pc (t) = dt {-i[H , p c {t)} + V[ Cl ]p c (t)} + dW{t)H[c x ]p c {t), 
where the nonlinear superoperator Ti is defined by 

TL[a]p = ap + pa) — Tr[(ap + pa')]p, 



(67) 



(68) 



and where dW is an infinitesimal real Weiner increment satisfying dW(t) 2 = dt. The subscript c indicates that the 
state is conditioned on the signal photocurrent, which is given by 



(69) 



where £,(t)dt = dW(t) represents Gaussian white noise |'2('|. These equations can be derived from the quantum jumps 
mentioned in Sec. Ill A, in the limit that the local oscillator amplitude used in the homodyne detection becomes 
infinitely large HQ]|. 

The feedback is then effected by introducing the time-dependent Hamiltonian 



H fb (t) = M c (t)Y. 



(70) 



Taking into account that the photocurrent in l|70[) is in reality smooth [unlike the mathematical expression (|69(l ]. and 
that its time argument must include a delay, one can derive the master equation (|66[) relatively straightforwardly. 



B. Langevin equation 



If the input to the source cavity is not a vacuum, then the technique of adiabatic elimination of Sec. IV A cannot 
be used. Unlike the intensity dependent feedback of Sec. Ill, it is nevertheless possible to obtain a sensible result 
for non-vacuum input using a Langevin equation approach. The difference arises because the quadrature flux of a 
broad-band squeezed or thermal state is a well-defined stochastic quantity, whereas the photon flux is not well defined 
(it tends to infinity in the limit of white noise). Because of the thermal and squeezed terms, it is necessary to use the 
explicit quantum equation (|2*3|) also including the time delay r. Adding the intercavity coupling gives the feedback 
Langevin equation for an arbitrary operator a 



da 



1 



(N + l)(2c\aci — ac\ci — c\cia) + N(2ciac\ — ac\c\ — cic\a) 



+M[cl[c\,a]] + M*[ Cl ,[c u al 



dt 



+1 



(N + l)(2cJ>ac2 — ac) 2 C2 — a^c^a) + N(2c2ac\ — ac2c\ — C2c\a) 



+M[4,[4,a]]+ M*[c2,[c2,a\] dt 
-\dB\d - dB lC \,a] - ^[dB\c 2 - dB 2 c\, a], +i[H + (c 2 + 4)3, a]dt, 



(71) 
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where the time argument of dB 2 — dB\ + c\dt is t — t, and all other operators are at time t. An operator a\ for the 
source cavity will obey the equation 



da 



1 



(N + \){2c\acx - ac\a - c\aa) + N(2dac\ 



ac\c\ 



c\c\a) 



dt 



+M[c[,[c[,a]]+M*[c u [c 1 ,a]\ 
-[dB\c\ — dBic\, a] + i[H + (c 2 + c\) J, a]dt. 



(72) 



Evidently, to obtain a Langevin equation for a source cavity operator involving no driven cavity operators, it is 
necessary only to adiabatically eliminate c 2 . From Eq. I|71l) . this obeys 



c 2 (t) = ~jc 2 (t) - ^jb 2 {t - r) - U(t). 
As before, if 7 3> 1, then there is no harm in replacing c 2 by the slaved value 



09 (*) = 



-2 

7f 



6 2 (t-r) 



iJ(t) 



V7 



(73) 



(74) 



providing that the resulting term in the Langevin equation is treated in the Stratonovich sense. That is to say, the 
effective feedback Hamiltonian 



H fb (t) = h(c 2 + c\)J = ih 



4(*-r)( 



2iJ{t) 



-2iJ{t) 



b 2 (t-r) 



(75) 



is to be treated in the same manner as the coupling Hamiltonian l|20() . The effective Hamiltonian (|75|l has the same 
form as Eq. (|2C)|> . but for the replacement of y / 72C2 by —iY [as defined above H65(l ], Thus it is possible to use the 
same procedure to derive a Langevin equation for a source cavity operator a in the limit t — > 0. From Eq. I|25|l . this 
is immediately seen to be 



da = (N + l) 
+N 
+M 



^(2c\a Cl - ocjcx - c| Cl a) - i[F, [Y, a}} + (i[Y, a] Cl + i^pT, a]) 

i(2ciac| - ac x c\ - cic\a) - [Y, a}} + (-i[Y, a]c\ - ic^Y, a}) 
1 r 



[cl[c\,a]} - [Y,[Y,a]] + 2i[c\,[Y,a]] 



1 



+M*-[[ci, [ci,o]] - [F, [!>]] - 2z[ Cl , [y, 1 



-[dB{(ci - iY) -dBi(c} +«y),a] +i[i/ ,a]dt, 

where is a true Ito increment. 

The master equation corresponding to this Langevin equation is now easy to derive. It is 

P = (N + vfacJp-iiY^^ + pcty+Nfactlp + iiY^lp + pc!}} 

+M [c[,p]] + i[Y, [alp]}} + M* |i[d, [d,p]] - i[y [d,p]] 

+ (2iV + 1 + M + M*)D[y]p - i[fr , p]- 



(76) 



(77) 



This is the general master equation for quadrature-based feedback in the presence of quantum white noise. For a 
vacuum input, it reduces to the master equation derived in the preceding section l|66|) . Recall that there are two 
derivations of Eq. (|66J) . from the all-optical model which we have just generalized for a non- vacuum input, and from 
an electro-optic model of feedback based on homodyne detection [y, ■ The correspondence between the two models 
is again close, as they can both be derived directly from a time dependent Hamiltonian. In the case of all-optical 
feedback, we have just seen that the effective feedback Hamiltonian (|75|l is 



H fb = h{b 2 + b\)Y. 



(78) 



This is identical to the current controlled Hamiltonian (|70|l if the homodyne signal is identified with the quadrature 
flux output of the cavity. 
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This close relationship suggests that it should be possible to derive the full master equation l|77|) from a measurement 
theory approach to homodyne measurement in the presence of noise. Such a theory has not been published before, 
to our knowledge. The difficulty in formulating the theory is that any white noise will, in theory, cause an infinite 
photon flux at the detectors. Thus, for any finite local oscillator strength, the signal will be swamped, just as it is 
for direct detection. A new approach is needed, in which the limit of infinite local oscillator strength is taken before 
determining the effect of the measurement on the system. After some effort, it can be shown that the generalization 
of the homodyne photocurrent expression (|69[1 is 

I(t)={ Cl + 4)(t) + VLt(t), (79) 

where we have defined 

L = 2N+1 + M + M*, (80) 

and £(t) is delta function normalized real Gaussian noise as before. Note that L can be less than one (its vacuum 
value) in the case of squeezed input noise. The stochastic evolution of the system given this photocurrent is 

{M M* 1 

(N+ l)V[ Cl ]p+ NV[c\]p+ — [cj, [c\,p\] + — [d, [ Cl ,p]}-i[H ,p]j 

dW(t) H [(N + M + l)c — (N + M *)c f ] p, (81) 

where the nonlinear superoperator TL is as defined in (|68|) . If the current (|79|l is fed back in the same manner as before 
[Eq. JTUJl], then a master equation can be derived using the same procedure as in the case of a vacuum input. The 
result is, not surprisingly, Eq. (|77|l . Thus, there is a complete correspondence between electro-optic and all-optical 
feedback in the case of quadrature-dependent feedback. 



C. In-loop and output squeezing 



As with intensity-dependent feedback, the all-optical quadrature-dependent feedback scheme described above is 
not a good way to try to produce nonclassical light in the source cavity. This is evident from the master equation 
(|64|l . Ignoring H , this will only produce squeezing if Y is a nonclassical operator (see preceding section). Given 
the origin of Y in the original coupling (|58|) . it is evident that at least a nonlinearity would be required to 
produce nonclassical source light. However, as noted in the preceding section, this requirement does not hold for the 
production of squeezing in the output light from the second cavity. Of course, this beam of light does not exist in 
the electro-optic case; it is only present for an all-optical feedback system. It turns out that this light can exhibit 
complete squeezing on-resonance, even though the source cavity is in a classical state. 

It is easiest to examine the in-loop and output fields in the Langevin equation approach. The output field operator 
&3 from the driven cavity is related to the in-loop field (denoted hi as above) by 

h = b 2 + ^c 2 . (82) 

Using the slaved value of c 2 J73J), this gives 

63 = - (6a - iY) = - (h + ci - iY) . (83) 

Apart from the unimportant sign change, this is just what would be expected from the effective Hamiltonian (|78|l . 

Now consider the simplest possible case, where b\ is in a vacuum state and Y is proportional to the imaginary 
quadrature of the source cavity 

Y = ~{-ic 1+ ic\). (84) 

This gives linear equations of motion by which the feedback alters the statistics of the real quadrature. Specifically, 
using x = ci + c\ and y = —ic\ + ic\, we get from Eq. I|76|l 

+ u x ^j - X(x + v x ) (85) 

1 

2 



~2/ + »y ) ■ (86) 
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Here v x and v y are vacuum quadrature noise operators which obey 



(Mt)Mt')) = Wy(t>v(t')) = S(t-t') (87) 



Since the equations for the two quadratures are uncoupled, the operator nature of v x and v y (as evidenced by Eq. 
is mostly unimportant. Hence it is easy to derive the steady state variance in x to be 

V{x)„ = <±±^. (89) 

This is always greater than one for any non-zero A, which is as expected since this sort of feedback cannot produce a 
squeezed intracavity state. 

To examine the extracavity fields, it is necessary to consider the noise at different frequencies. Defining the Fourier 
transform of a quantity by 

/oo 
e iut a(t)dt, (90) 
-oo 

Eqs. I|85I86[) become 

x(u) = - 1 ^" hA , v g (u) (91) 

y(cj) = - T l^v y {u). (92) 

The frequency domain counterparts to the time domain relationships (|87I88|) are identical but for the replacement of 
5(t — t') by 2ir5(ui + u r ). The in-loop quadratures of 62 are 

h + , n / \ 

i 2 (w)= = -TV^ (93 

\ + A — ILU 

fo(">) = -I^Ww). (94) 



2 " ZLJ 



These obey the commutation relations 



[£ 2 M, «2(w')] = 4 ^(^ + <->') x 1 ~ 2 i J _ + W 2 i T 7 ! ' ( 95 ) 

+ + A (i + «o;J 
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At first sight, the presence of the second factor here would seem to be a flaw in the theory, as it shows that the 
in-loop field does not obey the usual commutation relations for a free field. The commutator in Eq. I|95|) vanishes for 
low frequencies in the limit A — > 00. However, it must be remembered that the canonical commutation relations @ 
for the electromagnetic field are defined between the field at different points in space, but at the same time. If they are 
defined in terms of frequency components, as here, then these are strictly speaking spatial, not temporal, frequencies. 
It is only for a field which is free to propagate over an infinite space that the distinction vanishes. In our case, the 
in-loop field is spatially confined to a length cr, and we have let r go to 0. If we put r back in the equations, then it is 
not difficult to see that the only modification is to replace A in the frequency domain by Ae MT . With this replacement 
in the above frequency commutation relations, it is possible to show 4] that the time commutation relations are only 
changed for times greater than r. That is to say, the canonical (spatial) commutation relations are never violated. 
The field with operator & 2 (0 does not persist for longer than the time r as it travels from the first to the second 
cavity, and so there is never in existence two values of the field which violate the canonical commutation relations. 

Returning to the output field, we find from Eq. I|83|) that the quadratures of the output field 63 are (including the 
time delay r), 
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That is to say, the x quadrature of the output is unchanged from the 62 value (apart from the phase change due to 
the time delay), but the y quadrature has picked up an extra term. This extra term ensures that 

[x 3 (w),jk(<J)] = 4an6(u + u'), (98) 

as required because these fields may propagate to infinity. The spectrum for the real quadrature (which is equal to 
the spectrum of the photocurrent from a homodyne measurement) is defined by 



S x (u) = I (x(u)x(-ui'))duj', (99) 
and similarly for y. For the output field & 3 we get 



I ^ + Ae luJT — iw\ z 

SKM = i/s 3 ». (101) 

Note that for A — > 00, we get perfect squeezing at low frequencies for the x quadrature, and infinite noise in the y 
quadrature, as required by Heisenberg's uncertainty relations. 

Is there a simple way to understand this result, that the output may be perfectly squeezed, even though the cavity 
variance in x is unbounded [see Eq. l|89ft|? We look for an answer by an analogy with electro-optic feedback. Obviously 
the above analogy with electro-optic feedback mediated by homodyne detection will not suffice, because such detection 
destroys the output beam. What would be necessary would be a QND (quantum non-demolition) [2l| measurement 
of the output field. Feedback based on a QND measurement of output quadrature would have the same effect on 
the intracavity field as homodyne detection, provided it was efficient. It is not able to produce intracavity squeezing 
because an extracavity measurement is a poor measurement of the intracavity quadrature. However, if the output 
field x quadrature emerges from the QND device unchanged (or relatively little changed), then this quantity (the 
output field quadrature) can be well controlled by feedback. The fact that the feedback acts on the source cavity is 
relevant only so far as it affects the output. Such QND feedback schemes have been considered by Yamamoto et al 
0,0 and Shapiro et al 0. 

In the all-optical quadrature feedback scheme, the second cavity can be considered as a QND apparatus for the x 
quadrature output of the source cavity. As shown above, it does not alter the statistics of the x quadrature which 
reflects from it. However, it increases the variance of the output y quadrature, as required for a QND apparatus. The 
nonlinearity used in this case (coupling the quadrature of one mode to that of another) is precisely what has been used 
to model ideal QND quadrature measurements |?ll2l|. Rather than giving a current out as its measurement result, the 
output is directly coupled into the dynamics of the source cavity via the interaction Hamiltonian l|58|l . This enables 
fluctuations in the quantity it is monitoring (the output x quadrature) to be suppressed to an arbitrary degree. The 
Langevin equations for the all-optical case are identical to those of the electro-optic case, once the intervening current 
in the latter case has been eliminated. Thus we see that the correspondence between all-optical and electro-optic 
feedback can be made complete by using QND detectors, rather than normal (quantum demolition) photodetectors. 



V. COMPLEX AMPLITUDE FEEDBACK 



The two all-optical feedback schemes analyzed so far have both had an equivalent electro-optic scheme. The reason 
for this is that the couplings between driven and source cavity were QND couplings for the driven cavity. Since 
the driven cavity is slaved to the output of the source cavity, that means that the coupling effectively acts as a 
measurement of the output field of the source cavity, the result of which directly acts on the source cavity. In this 
section, we consider a Hamiltonian coupling between the two cavities which does not factorize as the direct product 
of a Hermitian operator in each cavity. We shall see that the feedback master equation obtained cannot be derived 
from any electro-optic scheme. 



A. General equations 



Consider the following Hamiltonian coupling between the two cavities 

V = h{c 2 B^ +4b), 



(102) 
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where B is an operator on the source cavity. If B is Hermitian (up to a phase factor), then this coupling is equivalent 
to that considered in the previous section on quadrature-dependent feedback. In general, however, this feedback is 
sensitive to both quadratures simultaneously and hence we have dubbed it complex amplitude feedback. Nevertheless, 
the analysis is identical to that used in the case of quadrature feedback. We will quote the main results. 

The master equation arising from complex amplitude feedback, including a squeezed or thermal input bath to the 
source cavity, is 



(N + 1) s V[ Cl +A]p-i 



2{c\A-A^ Cl ),p 



+M^[c\ + A\[c\+A\p]] + i 



- 1/ <i ^[ci + A, [d+A,p]] + i 



X+N^ct+A^p-i |( Cl A+-A4),p } 



A,A\p 



[d,A],p 



-i[H ,p] 



(103) 



where 



.4 



2iB 
7^' 



(104) 



where 7 is the large damping rate of the driven cavity as previously. It can be verified that in the case A = —iY, 
where Y is Hermitian, this equation is equal to Eq. (|77|l for quadrature feedback. 

One feature which distinguishes Eq. I|103|) from the quadrature feedback equation (|77|) is that it can produce a 
nonclassical state in the source cavity even if A is linear in c\ and c\. To see this, consider the case N = 0, Hq = 
to prevent any obscuring effects. The feedback master equation can be rewritten as 



p = D[ci]p + V[A] + (Apc\ + c lP A f - A^cip - pc[A). 



(105) 



For A linear, this equation can be converted into a Fokker-Planck equation for the Glauber-Sudarshan P function 
representation of the density operator . The condition for an initially positive P function (representing a classical 
state) to remain so is that the diffusion matrix be positive semi definite. The first term (damping) and third term 
(enclosed in round brackets above) will only give first order derivatives. Thus we need consider only the second term. 
Let 



A=-( Cl +pc{). 



(106) 



Then it can be readily shown that the eigenvalues of the diffusion matrix are proportional to \p\ 2 ± \p\. That is 
to say, if < < 1, then this complex amplitude feedback will produce a nonclassical state in the cavity. The 
= 1 limit gives the case of quadrature feedback, which, as we showed in Sec. IV C, cannot produce a nonclassical 
intracavity state. The other limit at \p\ — corresponds to the simple feedback considered in Sec. II B, which does 
not even require a nonlinear crystal (it is a classical geometrical optics problem). The property of nonclassicality 
distinguishes this all-optical feedback from any form of electro-optic feedback, because the latter is known not to 
produce nonclassical light from a linear feedback operator. 

The effective feedback operator (|106fl can be achieved from the interaction Hamiltonian 



V 



ihg (c 2 c\ 



c\d) 



p(c 2 d 



CnC 



■2°U 



(107) 



where we have put A and p real for simplicity, and where g — ^/t"A/4. Let the second mode C2 be an orthogonal 
polarization mode of the cavity containing the first mode d ■ That is to say, the second cavity is physically the same 
as the first, unlike the diagrammatic representation in Fig. 1. Then the first term in the Hamiltonian (|107f) could 
describe mode conversion, via a polarization rotator. The coupling constant g would be proportional to the (small) 
proportion of light converted at each pass, divided by the round trip time of the cavity. The second term, with 
strength gp, could only be produced by a nonlinear medium, such as a x' 2 ' crystal. The two polarization modes 
would be the signal and idler, and the second harmonic would have to be strongly driven and heavily damped so that 
it could be adiabatically eliminated. There are no obvious bars to setting up this scheme experimentally. 
The Heisenberg picture Ito equation equivalent to the master equation l|103|) is 



da 



(N + l) 



(^2c[ad 



ac\d 



c\da 



[c\A-A^d,a]) 



dt 
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•— (^2,c\ac\ — acic\ — Cic\a — [c\A^ — Ac\, a]^j dt 

[c{, [c\,a]} + [[clA^],a])dt+ — ([ci, [c u a] 
-\dB\c\ — dBic\, a] + i[Hq, a]dt, 



\d,A],a])dt 



where we have defined 



ci=ci+ A. 

This equation can be derived from the effective feedback Hamiltonian 

H fb = ih(b\A- A% 2 ). 
The output field, reflected off the mirror of the second cavity, is 

63 = -(&2 + -4) = -(&i + ci + A). 



(108) 
(109) 
(110) 
(111) 



Again, these equations are equivalent to the quadrature feedback equations when A = —iY. 

Consider the case where A is given by Eq. (|106fl . with A real and positive and [i real. This gives independent linear 
equations for the quadratures x, y of the field defined in Sec. IV C. Specifically from Eq. I|108J1 with N = M = 0, one 
obtains 



1 + A(l - /i) + 



(1-M 2 ) 



1 + ^(1-/4 



(112) 



where v x is as defined in Sec. IV C. The equation for the y quadrature is identical, but for the replacement of x by y, 
v x by v yi and fj, by — fj,. Note that for /j, = — 1, these equations agree with Eqs. (|85I86|I . as required. The intracavity 
steady state variance for x is 



V(x) 



l + A(l-/i) + (|) 2 (l-^)' 



(113) 



Note that for < /1 < 1, this implies a variance in x less than the unit variance of a coherent state. For < — n < 1, 
the variance in x will be greater than one, but that for y [obtained by replacing fj, by — fj, in Eq. (|113f) ] will be less 
than one. This is in accord with the result stated above, that for < \fi\ < 1, this linear all-optical feedback will 
produce a non-classical intracavity state. Furthermore, in the limit where A is very large, and e = 1 — fj, is very small 
(but not as small as A -1 ), then 



V{x) ss -> e/2. 



(114) 



That is to say, the intracavity state can be arbitrarily squeezed. This ideal result could preumably be obtained from 
electro-optic feedback, with some form of nonlinear feedback Hamiltonian. However the nonlinear nature of such 
feedback is quite different from the linear all-optical feedback considered here. 

Now consider the output squeezing. From the method of Sec. IV C, the output x spectrum can be calculated to be 



l(q + A M ) 2 +u, 2 



where 



1 + A 



2> t 1 -"")- 



(115) 



(116) 



The spectrum S^(u>) can be found by replacing /i by — /i, as before. As shown in Sec. IV C, the output of the feedback 
loop can show perfect squeezing even for 'classical' feedback with = 1. Consider < |^| < 1, so that a is always 
positive. For low frequencies, the output x quadrature is squeezed, with 5f (0) < 1, for fi < 0. Note that this is the 
sign of /1 which produces intracavity squeezing in the y quadrature. For /1 > 0, the output y quadrature is squeezed, 
while inside the cavity the x quadrature exhibits the nonclassical statistics. In understanding these counter-intuitive 
results, it must be remembered that the output beam 63 is not simply the output of the source cavity; the statistics 
of both quadratures are changed by its action as the feedback control beam. 
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B. Electro-optic analog 



As stated above, all-optical complex amplitude feedback has no electro-optical counterpart in general. This is 
because it is not possible to measure both the real and imaginary quadratures of the output field simultaneously with 
unit efficiency. Even a QND measurement of one quadrature would introduce noise into the other and so prevent a 
measurement of both. However, it is possible to do two inefficient measurements of both quadrature with the two 
efficiencies adding to one (or less than one in practice). It is simplest to consider heterodyne detection, which is 
equivalent to a homodyne measurement of each quadrature, each with efficiency of one half. The effect of inefficient 
measurement is to increase the noise introduced by the feedback by a factor inversely proportional to the efficiency. 
The precise meaning of this statement will become clear shortly. 

In order to compare this all-optical complex amplitude feedback to feedback from heterodyne detection, it is 
convenient to rewrite Eq. (1 103|) in terms of the Hermitian operators X, Y defined by 

A = X — iY. (117) 

We obtain 

p = -i[H , P ] + {N + 1) \v[ Cl ]p - i[Y, c lP + pc[} - i[X, -ic lP + ipc\] + V[X - iY]pj 
+N \v[c\]p + i[Y, c\p + pci] - i[X, ic\p - ipa] + V[X + iY]p\ 

+M {-[4, [c\, p]} + i[Y, [c\, p\] - i[X, [ic\, p\] -\[X + iY, [X + iY, p}] 

+M* |i[ci, [ci,p]] - i[Y, [ Cl ,p]] + i[X, [-id, p]] + \[X- iY, [X - iY, p]]| . (118) 

If A is linear in the field amplitude, then X and Y are proportional to two orthogonal quadratures of the field. For 
controlling noise, it would be sensible for these to be proportional to the x and y quadratures respectively. 

The terms in this equation linear in the feedback operators X and Y are reminiscent of the terms in the quadrature 
dependent feedback master equation \77f. However, the terms bilinear in X, Y, which we will refer to as noise terms 
for reasons which will become evident later, are different. As noted in the preceding section, quadrature feedback can 
be effected by either an all-optical scheme or a homodyne detection scheme. In the latter case, the master equation 
was derived from the stochastic evolution equation describing the effect of homodyne detection on the system. The 
equivalent equation for heterodyne detection is found to be 



dp(t) 



( M M* 1 

= dtUN + l)V[c 1 ]p + NV[c\]p + — [c\, [c\,p}\ + — [a, [d,p]] -i[H 0)P n 

+—^dW x (t)H [(N + M + l)c-(N + M*)J] p 

+ ^^dW y (t)H [(N — M + l)(—ic) - (N-M*)(iJ)] p, (119) 



where 

L x = 2N+1 + M + M* (120) 
L y = 2N+1-M-M*. (121) 

Here, dW x (t) and dW y (t) are real infinitesimal Weiner increments. Physically, they are related to the noise in 
the two photocurrents which arise from the heterodyne detection by dW q (t) = t; q (t)dt (q = x, y). These currents are 
given by 

/*(*) = (ci+c{) + ^2L x ~t x (t) (122) 
I y {t) = (-id +ic\) + y/2L^S y (t), (123) 

where the normalization has been chosen to make the deterministic term the same as in the homodyne case. Note 
that the noise is greater in this case, with L q multiplied by 2. This is because each measurement is of efficiency half. 
In general, if the efficiency of the extra-cavity measurement of the quadrature q are r] q (such that rj x + rj y < 1), then 
2L q is replaced by L q jr\ q in the above equations. 
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Now the stochastic equation 1)119)1 can be used to derive a feedback master equation analogous to the all-optical 
Eq. 1)11811 . The currents control the time-dependent feedback Hamiltonian 

H fb (t) = h[I x (t)Y + I y (t)X}. (124) 

This should be compared with Eq. 1)11(1 f) . It is not equivalent to that equation, because there is more noise in the 
currents than there is in the quadratures. If one wished to represent the currents by operators as was done for the 
case of homodyne detection, then this extra noise would enter at the beam splitter necessary to split the output into 
two separate homodyne devices. By adding the evolution from the Hamiltonian (1 124)1 to that of Eq. ()119)l . we obtain 
the heterodyne feedback master equation 

P = {N + l){v[a]p- i[Y,c lP + pc\]- i[X,-icip + ipc\]} 
+N VD[c\]p + i[Y, c\p + pci] - i[X, ic\p - ipd]\ 

+M ji[cj , [c\,p]] + i[Y, [c\, p}} - i[X, [ic\,p]] 

+M* |i[ci, [<n,p]] - i[Y, [ci, p]] + i[X, [-icL,p]] 

+2L x V[Y]p + 2L y V[X]p - i[H , p]. (125) 

Note that the desired feedback terms (linear in X and Y) are the same in this equation as in the all-optical Eq. ()118)l . 
but the diffusion terms are different. 

To elucidate this difference, we return to the most basic example of all-optical feedback, considered in Sec. II B. 
The output of the source cavity is simply fed back into another mirror of that cavity. Such feedback is covered by the 
master equation derived in this section, with 

A = yfyeta. (126) 

Here, a is the annihilation operator for the source cavity with decay rate 7. In this simple case, there is no need for 
a second cavity. As derived in Sec. II, the master equation for the cavity is 

p = 27(1 + cos cj))V[a]p - ij sin (j)[a^ a, p), (127) 

where the input is in the vacuum state. Attempting to replicate this feedback by using heterodyne detection yields 
the master equation 

p = 27(1 + cos (f>)V[a\p - i7sin</>[a l 'a, p]+ 1 (V[a]+V[a)])p. (128) 

The equation of motion for the mean field from this equation is identical to that of Eq. 1)127)1 . However, the presence 
of the extra term introduces noise into both quadratures equally. If 4> = n, so that the deterministic dynamics are 
eliminated, then the variance in each quadrature will simply grow linearly. This clearly shows the effect of the noise 
introduced by attempting to measure both quadratures in electro-optic feedback, as opposed to the coherent back 
coupling of both quadratures in all-optical feedback. 



VI. CONCLUSION 



At the simplest level, feedback is a process by which a system influences itself through its action on a second system. 
Usually, the second system would be a complicated feedback apparatus, consisting of measurement devices, signal 
processors and the like. However, it is always possible to conceive of more direct schemes, in which the feedback loop 
is treated on the same level as the system. In this paper we have examined feedback on optical cavities. The usual 
implementation of optical feedback is electro-optic feedback. The light emitted by the cavity enters a detector, and 
the photocurrent produced is used to control the dynamics of the cavity by some electro-optic devices. The more 
direct method is all-optical feedback. We have considered turning the output beam from the source cavity onto a 
second cavity which is directly coupled to the first by some nonlinear crystal. To compare these two methods at the 
quantum limit, we assumed that the feedback could be approximated as a Markovian process. In the electro-optic 
case, this corresponds to assuming that the overall time delay in the feedback loop is much smaller than the lifetime 
of the source cavity, while in the all-optical case, we also need the linewidth of the second cavity to be much greater 
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than that of the first. Under these conditions, a master equation, or quantum Langevin equation can be derived for 
the source cavity alone. 

Our main conclusion is that there is a strong relationship between all optical and electro-optical quantum-limited 
feedback. In fact, if the direct coupling Hamiltonian V between the two cavities factorizes as the product of Hcrmitian 
operators in each cavity, then the all optical feedback has an exact electro-optic feedback counterpart, at least as far 
as the effect on the source cavity is concerned. If V depends on the intensity of the driven cavity, then the result 
can be reproduced by feedback based on direct photodetection, while quadrature-dependent all-optical feedback can 
be reproduced by homodyne detection. The explanation for this is that the nonlinear interaction acts as an effective 
measurement of the driven cavity intensity or quadrature, which (because it is heavily damped) is a measurement of 
the intensity or quadrature of the output of the first cavity. Unlike the electro-optic case, the measuring device then 
directly influences the source cavity through the same coupling, instead of having all of the intermediate equipment. 
On the other hand, a purely optical feedback scheme which is sensitive to both quadratures of the driven cavity 
cannot be replicated by currents. That is because such an interaction is not like a simultaneous measurement of both 
quadratures of the driven cavity. Rather, the complex amplitude of the driven cavity interacts coherently with another 
non-Hermitian quantity in the source cavity. An analogous electro-optic scheme using heterodyne detection can be 
defined, which has the same semiclassical effect as the complex amplitude feedback. However, the two independent 
measurements of the two quadratures make the feedback incoherent and introduces extra diffusion terms into the 
master equation. 

The equivalence of electro-optic and all-optical feedback for the intensity and quadrature feedback means that the 
latter is subject to the same restrictions as the former as far as the generation of nonclassical light is concerned. 
Specifically, if the nonlinear interaction Hamiltonian V is linear in the source cavity amplitude or intensity, then the 
feedback cannot produce a nonclassical state in the source cavity. The electro-optic equivalent of this theorem is 
that feedback based on extra-cavity detection cannot produce nonclassical light by driving or detuning the cavity. 
However, there is a property of all-optical feedback not present in electro-optical schemes which makes the situation 
less clear cut: an output beam. This is the feedback loop beam reflected from the second cavity. It turns out that 
this beam may be arbitrarily squeezed, even though the source cavity remains classical. 

The way to understand this is using the measurement analogy explained above. The driven cavity is a QND 
measurement device for say the x quadrature of the output of the source cavity. This device also may control the 
dynamics of the source cavity such that the measured fluctuations in the x quadrature are suppressed. Thus the x 
quadrature of the output of the source cavity can be squeezed; its statistics are unchanged upon reflection at the 
second cavity. However, it is not until the output reflects off the second cavity (the QND device) that the extra 
fluctuations in the y quadrature are put in, due to the effect of the measurement. This would seem to indicate that in 
the loop, Hciscnbcrg's uncertainty relations fail, even though they are satisfied in the output loop (which is all that is 
observable). However, this is not the case. The canonical commutation relations, which are actually between different 
spatial points of the field at the same time, are never violated. The anticorrclations in the output field which result 
in squeezing are only present between parts of the field which are separated in time by more than the time delay in 
the feedback loop. There is no contradiction. 

Thus, it would seem that there is one potential application for all-optical feedback: producing squeezed light. 
However, in order to build such a device producing well-squeezed light, the coupling would have to use a frequency- 
(but not polarization-) degenerate \^ crystal, as well as a polarization converter. It would seem easier to use 
the traditional squeezer, a x crystal acting as a degenerate parametric oscillator. Intensity-dependent all-optical 
feedback is even less practical, requiring a low loss x^ nonlinearity to operate. The smallness of higher order 
nonlinearities is sufficient justification as to why we have not considered all-optical feedback with a coupling dependent 
on higher order field moments of the driven cavity. In fact, such higher order feedback does not produce any new 
results. At least in the regime where the second cavity can be adiabatically eliminated, the higher order terms either 
give a vanishing contribution, or reproduce the results of amplitude or intensity feedback. Thus we can conclude 
that all-optical feedback is probably not a practical way of controlling quantum noise, although there may be other 
applications. Nevertheless, the predicted results are interesting, and some experiments should be feasible with current 
technology. The similarities with and differences from electro-optic feedback yield important insights about the nature 
of feedback in general. 
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